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Abstract 

In this paper, we present topological black holes of third order Lovelock gravity in the presence 
of cosmological constant and nonlinear electromagnetic Born-Infeld field. Depending on the met- 
ric parameters, these solutions may be interpreted as black hole solutions with inner and outer 
event horizons, an extreme black hole or naked singularity. We investigate the thermodynamics of 
asymptotically flat solutions and show that the thermodynamic and conserved quantities of these 
black holes satisfy the first law of thermodynamic. We also endow the Ricci flat solutions with a 
global rotation and calculate the finite action and conserved quantities of these class of solutions 
by using the counterterm method. We compute the entropy through the use of the Gibbs-Duhem 
relation and find that the entropy obeys the area law. We obtain a Smarr-type formula for the 
mass as a function of the entropy, the angular momenta, and the charge, and compute temperature, 
angular velocities, and electric potential and show that these thermodynamic quantities coincide 
with their values which are computed through the use of geometry. Finally, we perform a stability 
analysis for this class of solutions in both the canonical and the grand-canonical ensemble and show 
that the presence of a nonlinear electromagnetic field and higher curvature terms has no effect on 
the stability of the black branes, and they are stable in the whole phase space. 
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I. INTRODUCTION 



Over the last few years, several extra-dimensional models have been introduced in an 
attempt to deal with the hierarchy problem. These models can lead to rather unique and 
spectacular signatures at Terascale colliders such as the LHC and ILC In higher dimensions, 
it is known that the Einstein-Hilbert (EH) Lagrangian, R, can only be regarded as the first 
order term in an effective action, so one may on general grounds expect that as one probes 
energies approaching the fundamental scale, significant deviations from EH expectations are 
likely to appear. This motivates one to consider the more general class of gravitational 
action: 

/g= lib/ dn+lx ^9HR,R^ u ,R^ P ^ pa )- 

The presence of higher curvature terms can also be seen in the renormalization of quantum 
field theory in curved spacetime [l| , or in the construction of low energy effective action of 
string theory [2] . Among the higher curvature gravity theories, the so-called Lovelock gravity 
is quite special, whose Lagrangian consist of the dimensionally extended Euler densities. 
This Lagrangian is obtained by Lovelock as he tried to calculate the most general tensor 
that satisfies properties of Einstein's tensor in higher dimensions Since the Lovelock 
tensor contains derivatives of metrics of order not higher than two, the quantization of 
linearized Lovelock theory is free of ghosts 4|. Thus, it is natural to study the effects of 
higher curvature terms on the properties and thermodynamics of black holes. 

Accepting the nonlinear terms of the invariants constructed by Riemann tensor on the 
gravity side of the action, it seems natural to add the nonlinear terms in the matter action 
too. Thus, in the presence of an electromagnetic field, it is worthwhile to apply the action 
of Born-Infeld |J instead of the Maxwell action. In this paper, we generalize static and 
rotatingblack hole solutions of third order Lovelock gravity in the presence of Maxwell 
field |6|, |7| to the case of these solutions in the presence of nonlinear electromagnetic fields. 
Indeed, it is interesting to explore new black hole solutions in higher curvature gravity 
and investigate which properties of black holes are peculiar to Einstein gravity, and which 
are robust features of all generally covariant theories of gravity. The first aim to relate 
the nonlinear electrodynamics and gravity has been done by Hoffmann [§]. He obtained a 
solution of the Einstein equations for a pointlike Born-Infeld charge, which is devoid of the 
divergence of the metric at the origin that characterizes the Reissner-Nordstrom solution. 
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However, a conical singularity remained there, as it was later objected by Einstein and 



Rosen. The spherically symmetric solutions in Einstein-Born-Infeld gravity with or without 
a cosmological constant have been considered by many authors 9I. llo|. while the rotating 
solutions of this theory is investigated in Also, these kinds of solutions in the presence 
of a dilaton field have been introduced in 12j. The static black hole solutions of Gauss- 
Bonnet-Born- Infeld gravity have been constructed in Ref. [3], and the rotating solution of 
this theory has been considered in }14| . 

The out line of our paper is as follows. We present the topological black holes of third 
order Lovelock gravity in the presence of Born- Infeld field in Sec. [Ill In Sec. IHIl we calculate 
the thermodynamic quantities of asymptotically flat solutions and investigate the first law 
of thermodynamics. In Sec. IIVI we introduce the rotating solutions with flat horizon and 
compute the thermodynamic and conserved quantities of them. We also perform a stability 
analysis of the solutions both in canonical and grand canonical ensemble. We finish our 
paper with some concluding remarks. 

II. TOPOLOGICAL BLACK HOLES 

The action of third order Lovelock gravity in the presence of nonlinear Born-Infeld elec- 
tromagnetic field is 

Ig = l^r" / dn+lx ^9 (-2A + A + oc 2 C 2 + a 3 C 3 + L(F)) , (1) 

where A is the cosmological constant, a 2 and a 3 are the second and third order Lovelock 
coefficients, £1 = R is just the Einstein-Hilbert Lagrangian, C 2 = Rfi^sR 11 ^ 6 — ^R^R^+R 2 
is the Gauss-Bonnet Lagrangian, 

C O TipVOK T3 JDPT I Q T}P V T3 aK TDpT | r) A TJ/IVCTK p DP 

+3RR^ K R^ U + 2AR^ K R^R KU + WR^R ua R\ - YHUTR^ + R 3 (2) 
is the third order Lovelock Lagrangian, and L(F) is the Born-Infeld Lagrangian given as 



JW=v(l-^l + £). (3) 

In the limit (3 — > 00, L(F) reduces to the standard Maxwell form L(F) = —F 2 , where 
Ffiu — d^A u — d u A fM . Varying the action (frj with respect to the metric tensor g^ u and 



electromagnetic vector field Ap the equations of gravitation and electromagnetic fields are 
obtained as: 

' '' (4) 



Gg) + Ag, u + a 2 G$ + a 3 G$ = \g, v L{F) + 



1 + H 

1 + 2/3 2 



<9„ 




0. 



(5) 



where G$ is the Einstein tensor, and Gjw and are the second and third order Lovelock 
tensors given as (la ]: 

= ^{R^ktRJ 7 ^ — ^R^puaR pCT — 2R[i<jR a v + RR^u) — -£29^, (6) 



/t(3) O/ /I JDTpCFK p pA QD T /' E> CTK t?A inn TOKT3 TD^P 

(J >f ~~ o^n, n aK \ p n UTfl on Acr n r/1 .n -r ^-f^, JX^KXpn, T ^ 

_ TiTpCfK p D _l_ Q pT F?/ 3 _|_ S E" 7 F? T '° f? K 

11 ■ r hrKTp r h/lJ, 1 orL uap JX T/i- rt K ' 0rL urn 171 ap 1 ^ p 
1 A p ran p pp t<tk p pp . J prpcrK p p . 9 p p Krp p 

+8R T vppR p aR a T - %R C 7 ur P R' uRp - 8R T apR c ' T Ry P - ^RR T V ppR p T 

+AR Tp R pr R U p — 8R T u R Tp R p + ARR up R p — i? i?^) — -C 3 gp U . 



(7) 



Here we want to obtain the (n + l)-dimensional static solutions of Eqs. (jl]) and (jSJ). We 
assume that the metric has the following form: 

dr 2 



where 



dn 2 



ds 2 = -f(r)dt 2 + — — + r 2 dtt 2 , 
f( r ) 



del + E II sin2 d i de i 



(8) 



k = 1 



i=2 j=l 

d6f + s in h 2 e x de\ + sinh 2 B x E fi sin2 e j de I k = -1 

i=3 j=2 

ra-1 

E = o 

i=i 

represents the line element of an (n — l)-dimensional hypersurface with constant curvature 
(n — l)(n — 2)fc and volume Ki-i- 

Using Eq. (jSJ), one can show that the vector potential can be written as 



(9) 
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where q is an integration constant which is related to the charge parameter and 

_{n- \){n-2)q 2 
^ ~ 2f3 2 r 2n ~ 2 ' 

In Eq. (jHJ) and throughout the paper, we use the following abbreviation for the hypergeo- 
metric function 



1 n- 2 



3n-4 



_2'2n-2j ' [2n-2 
The hypergeometric function F(t])— ► 1 as rj — > (/? 

reduces to the gauge potential of Maxwell field. One may show that the metric function 



= F(z). (10) 
oo) and therefore of Eq. 



^-(l-^r) 1 / 3 ), 



(?(r) = 1 + 



3am 



'1 



satisfies the field equations 



f(r) = 

12a[3 2 
n{n — 1) 
in the special case 
a 

a 2 



A 



(n — 1)?7 
2/3^ (n-2) 



F(v) 



(11) 
(12) 



(n-2)(n-3)' 
a 2 



a 3 



72( n - A 2 y 



where m is the mass parameter. Solutions of Gauss-Bonnet gravity are not real in the whole 
range < r < 00 and one needs a transformation to make them real [l^. But, here the 
metric function f(r) is real in the whole range < r < 00. 

In order to consider the asymptotic behavior of the solution, we put m = q = where 
the metric function reduces to 



6Aa 



1/3- 



a \ n{n — I) J 

Equation (fT3"|) shows that the asymptotic behavior of the solution is AdS or dS provided 
A < or A > 0. The case of asymptotic flat solutions (A = 0) is permitted only for k — 1. 



As in the case of black holes of Gauss-Bonnet-Born-Infeld gravity [13|, |14| . the above 
metric given by Eqs. ([H]), ( TTTT) and ( |T2l) has an essential timelike singularity at r = 0. 
Seeking possible black hole solutions, we turn to looking for the existence of horizons. The 
event horizon(s), if there exists any, is (are) located at the root(s) of g rr = f(r) = 0. 
Denoting the largest real root of f(r) by r + , we consider first the case that /(r) has only 
one real root. In this case f(r) is minimum at r + and therefore /'( r +) = 0- That is, 

l)k [3(n - 2)4 + 3(n - 4)kar\ + (n - 6)k 2 a 2 ] + \2r\0 2 (l - y/l + 77+) - 6Ar 6 + = 0. 

(14) 



n 



One can find the extremal value of mass, m ex t, in terms of parameters of metric function 
by finding r + from Eq. (j!4p and inserting it into equation f(r+) = 0. Then, the metric of 
Eqs. (151) , (II ip and (TT2"j) presents a black hole solution with inner and outer event horizons 
provided m > m cxt , an extreme black hole for m = m ext [temperature is zero since it is 
proportional to /'(r + )] and a naked singularity otherwise. It is a matter of calculation to 
show that m cxt for k = becomes 

_ 2(n - l)g c 2 xt / A(A-4/3 2 ) \ (- 2 )/( 2 - 2 ) A(A-4/3 2 ) 
CXt " n ^2(n-l)(n-2)/3^ xt y' 1 40* j ' 

The Hawking temperature of the black holes can be easily obtained by requiring the ab- 
sence of conical singularity at the horizon in the Euclidean sector of the black hole solutions. 
One obtains 

/' (r+) in - l)k [3(n - 2)r| + 3(n - A)kar\ + (n - 6)k 2 a 2 ] + 12r 6 + /3 2 (1 - y/l + rj+) - 6Ar 



4tt 127r(n - l)r + (r 2 + kaf 

(16) 



It is worthwhile to note that T + is zero for m = m ext . 

III. THERMODYNAMICS OF ASYMPTOTICALLY FLAT BLACK HOLES FOR 

k = 1 

In this section, we consider the thermodynamics of spherically symmetric black holes 
which are asymptotically flat. This is due to the fact that only the entropy of asymptotically 
black holes of Lovelock gravity is well known [171 ] . Usually entropy of black holes satisfies 
the so-called area law of entropy which states that the black hole entropy equals one-quarter 
of the horizon area fl8| . One of the surprising and impressive features of this area law of 
entropy is its universality. It applies to all kinds of black holes and black strings of Einstein 



gravity 



in general 



^However, in higher derivative gravity the area law of entropy is not S ati sfl ed 
201 ] . It is known that the entropy of asymptotically flat black holes of Lovelock 



gravity is [17 ] 



[(n-2)/2] 



S=~ J2 kak S d n - l x^gC k -i, (17) 
k=i J 

where the integration is done on the (n — l)-dimensional spacelike hypersurface of the Killing 
horizon, g^ u is the induced metric on it, g is the determinant of g^, and Ck is the fcth order 
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Lovelock Lagrangian of g^ v . Thus, the entropy for asymptotically flat black holes in third 
order Lovelock gravity is 

S = \j d n ~ l x^ (l + 2a 2 R + 3a 3 (R^ aK R^ K - m^BT + R 2 )) , (18) 

where R^ U p a and R^ u are Riemann and Ricci tensors and R is the Ricci scalar for the induced 
metric g a b on the (n — l)-dimensional horizon. It is a matter of calculation to show that the 
entropy of black holes is 

Vn-i (a . 2(n-l) o , n-1 2 \ 5 



4\ + n - 3 + n - 5 / + 
The charge of the black hole can be found by calculating the flux of the electric field at 
infinity, yielding 

Q ^^ES, (20) 

The electric potential $, measured at infinity with respect to the horizon, is defined by 

$ = I™ - A,x% =T+ , (21) 
where \ = d/dt is the null generator of the horizon. One finds 



*^/|^f^''+>- < 22) 

The ADM (Arnowitt-Deser-Misner) mass of black hole can be obtained by using the behavior 
of the metric at large r. It is easy to show that the mass of the black hole is 

M=-n=±( n -i)m. (23) 

We now investigate the first law of thermodynamics. Using the expression for the entropy, 
the charge, and the mass given in Eqs. ( fT9l) . ( l20l) and ( |23|) . and the fact that f(r + ) = 0, one 
obtains 



M(S,Q) 



(n-1) 



2r\ 



16tt n(n — 1, 



2{3 2 



l- + 3 + 



in 



in 



1)9 



r(3) 



A 



r r r 2 + arT* - 



n— 4 



OCT 



+ 



(24) 



where 



16n 2 Q 2 
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2„2n-2 • 



/3 2 r 
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In Eq. (124)) r + is the real root of Eq. (fl9l) which is a function of S. One may then regard the 
parameters S and Q as a complete set of extensive parameters for the mass M(S, Q) and 
define the intensive parameters conjugate to them. These quantities are the temperature 
and the electric potential 

'd_M\ (d_M\ 
dS) Q \dQj s - 



T 



(25) 



Computing dM/dr + and dS/dr + and using the chain rule, it is easy to show that the 
intensive quantities calculated by Eq. ( 1251) coincide with Eqs. ( fl~6j) and ( 1221) . respectively. 
Thus, the thermodynamic quantities calculated in Eqs. (fl6l) and (1221) satisfy the first law of 
thermodynamics, 

dM = TdS + $dQ. (26) 



IV. THERMODYNAMICS OF ASYMPTOTICALLY ADS ROTATING BLACK 
BRANES WITH FLAT HORIZON 



Now, we want to endow our spacetime solution ([8]) for k = with a global rotation. 
These kinds of rotating solutions in Einstein gravity have been introduced in 21J. In order 
to add angular momentum to the spacetime, we perform the following rotation boost in the 
t — (pi planes 

t^->Et - cii<f)i, (pi h-> E(pi - -^t (27) 

for i = l...[n/2], where [x] is the integer part of x. The maximum number of rotation 
parameters is due to the fact that the rotation group in n + 1 dimensions is SO(n) and 
therefore the number of independent rotation parameters is [n/2]. Thus the metric of an 
asymptotically AdS rotating solution with p < [n/2] rotation parameters for flat horizon 
can be written as 



ds 2 



-f(r) Edt - j> 



+ 



dr 2 

W) 



8=1 
2 P 



2 P 

+ p-> , (aidt - El 2 d(pi)' 

8=1 

n-1 

— y^ y (aid<f>j ~ ajd(pi) 2 + r 2 dfc, 

i<j i=p+l 



(28) 



where E = 
written as 



1 + E! ia 2 /l 2 . Using Eq. ©, one can show that the vector potential can be 

(29) 



4* = - 



n 



1) ? 



2n — 4 r 



n-1 



F(v) - Slai) ( 



no sum on i) 



One can obtain the temperature and angular momentum of the event horizon by analytic 
continuation of the metric. One obtains 

f(r+) 



4vrS 2(n 



±— (2(3 2 (1 - v^T^) - A) , (30) 

a = ^ ( 31 ) 
where 7] + = rj(r = r + ). Next, we calculate the electric charge and potential of the solutions. 
The electric charge per unit volume V n -i can be found by calculating the flux of the electric 
field at infinity, yielding 

Using Eq. (j2ip and the fact that % = d t + Yli is the null generator of the horizon, the 
electric potential $ is obtained as 



A. Conserved quantities of the solutions 

Here, we calculate the action and conserved quantities of the black brane solutions. In 
general the action and conserved quantities of the spacetime are divergent when evaluated on 
the solutions. A systematic method of dealing with this divergence for asymptotically AdS 
solutions of Einstein gravity is through the use of the counterterms method inspired by the 



anti-de Sitter conformal field theory (AdS/CFT) correspondence 22] . For asymptotical!} 



AdS solutions of Lovelock gravity with flat boundary, R a bcd(l) = 0, the finite action is 



I = I G + ^~ I d n x V=7 {Lu + a 2 L 2b + a 3 L 3b } + — I d n x^f^ ( ^— - j , (34) 
8vr JdM 87r JdM V L J 



where L is 



15iya(l-A) 
5/ 2 + 9a-/ 2 A 2 -4/ 2 A' 1 ' 

A = (1 - ^f>\ (36) 

One may note that L reduces to I as a goes to zero. The first integral in Eq. f)34p is a 
boundary term which is chosen such that the variational principle is well defined. In this 
integral L lb = K, L 2b = 2( J - 2G { ^ K ab ) and 

L 3b = 3(P - 2G^K ab - 12R ab J ab + 2RJ - AKR abcd K ac K bd - 8R abcd K ac K b K ed ), 



where 7^ and K are induced metric and trace of extrinsic curvature of boundary, G^j and 
G^ b ' are the n- dimensional Einstein and second order Lovelock tensors (Eq. ([6])) of the 
metric 7 a b and J and P are the trace of 

J- = \(2KK^Kl + K c!t K cd K ab - 2K^K M - (37) 

and 

P afe = i{[if 4 -6ir 2 ir d K cd + 8^ c ^^ 

-(AK 3 - 12iTK ed K ed + 8K de K e f K fd )K ac K c b - 2AKK ac K cd K de K e b 

+ (12K 2 - 12K ef K £ f)K ac K cd K db + 2AK ac K cd K de K e f K bf }. (38) 

Using Eqs. (CD) and (IMjl . the finite action per unit volume V n -\ can be calculated as 

T + \l6nP 4n(n-l)7r 8n7rrf~ 2) J ' W 

I — I 

Using the Brown- York method [24]) the finite energy-momentum tensor is 

1 77 — 1 

T ab = _{( K ab _ Rl ab^ + 2a2 (3J ab - J 7 ab ) + 3a 3 (5P ab - P 7 ab ) + 7 af> }, (40) 

87T L 

and the conserved quantities associated with the Killing vectors d/dt and d/d(f) % are 

M = —m (nZ 2 - l) , (41) 
16tt V ; ' V ; 

J; = nEmai, (42) 

which are the mass and angular momentum of the solution. 
Now using Gibbs-Duhem relation 

k 

( M — Q<b — 1 

T 



1 k 
S = — (M - Q$ - ^fiiJi) - I, (43) 



and Eqs. ([33]), ([39]) and gTJ-gSJ one obtains 

S = |r»" 1 (44) 

for the entropy per unit volume V^_i. This shows that the entropy obeys the area law for 
our case where the horizon is flat. 
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B. Stability of the solutions 



Calculating all the thermodynamic and conserved quantities of the black brane solutions, 
we now check the first law of thermodynamics for our solutions with flat horizon. We obtain 
the mass as a function of the extensive quantities S, J, and Q. Using the expression for 
charge mass, angular momenta and entropy given in Eqs. (|32|) . (l4Tj) . (1421) . (1441) and the fact 
that f(r + ) = 0, one can obtain a Smarr-type formula as 



where J = |J| = yY^i J? an d Z = E 2 is the positive real root of the following equation: 

Z l/2(n-l) ^ [4/2/3 2 (1 _ V /T + ^) + n (n - 1)] gn/(n-l) ^^(^g) 

VZ^T ~ (n - l)ulJ + (n - 2) J(45)(™- 2 )/("- 1 ) ' (46) 

One may then regard the parameters S, Jj's, and Q as a complete set of extensive parameters 

for the mass M(S, J, Q) and define the intensive parameters conjugate to them. These 
quantities are the temperature, the angular velocities, and the electric potential 

dM\ fdM\ ^ (dM\ 

Straightforward calculations show that the intensive quantities calculated by Eq. (|47p coin- 
cide with Eqs. fl30l) . (13T1) and fl33|) . Thus, these quantities satisfy the first law of thermody- 
namics: 

k 

dM = TdS + J^JWi + $<2Q. 

i=i 

Finally, we investigate the local stability of charged rotating black brane solutions of 
third order Lovelock gravity in the presence of nonlinear electrodynamic Born-Infeld field in 
the canonical and grand canonical ensembles. In the canonical ensemble, the positivity of 
the heat capacity Cj ; q = T + /(d 2 M/dS 2 )j,Q and therefore the positivity of (d 2 M/dS 2 )j t Q is 
sufficient to ensure the local stability. Using the fact that 



3 3n-4 



5n — 6 
2n - 2 



n-l)z { y (3 2 S 2J y/TT 



2' 2n-2 
it is easy to show that 

d 2 M _ 2[(n - l)(n - 2) 2 q 2 - Arg^Vl±j+ + 2/3 2 r { + n ~ 2) + 1)] 
~dS 2 ~~ (n-l^rf^VTT^: 

8(~ 2 - 1)A (-r 2 "- 2 v ^T^(^ ii + 2/5 2 / 2 ) + (n - l)(n - 2)/5/ 2 g 2 + 2/3 3 Z 2 rf 
7rm/ 2 ^ 2 S 2 (n - 1) 4 (4S 2 + 1)(1 + r/+)r( 4 ™- 6 ) 
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{48) 



Both of the two terms of Eq. (14"8]) are positive, and therefore the condition for thermal 
equilibrium in the canonical ensemble is satisfied. 

In the grand canonical ensemble, the positivity of the determinant of the Hessian matrix of 
M(S, Q, J) with respect to its extensive variables X i7 H^. x . = (d 2 M/dXidXj), is sufficient 
to ensure the local stability. It is a matter of calculation to show that the determinant of 

a S,Q,J 1S - 



[n - 2){n - \) 6 mVc?r+ ' [(n - 2)^ 2 + 1\ ^1 + r] + r + 

(49) 

Equation (1491) shows that the determinant of the Hessian matrix is positive, and therefore 
the solution is stable in the grand canonical ensemble too. The stability analysis given 
here shows that the higher curvature and nonlinear Maxwell terms in the action have no 
effect on the stability of black holes with flat horizon, and these kinds of black holes are 
Lermodynamically stable as in the case of toroidal black holes of Einstein-Maxwell gravity 



251 ] . This phase behavior is also commensurate with the fact that there is no Hawking-Page 



transition for a black object whose horizon is diffeomorphic to W and therefore the system 
is always in the high temperature phase [26] . 



V. CLOSING REMARKS 



In this paper we considered both the nonlinear scalar terms constructing from the cur- 
vature tensor and electromagnetic field tensor in gravitational action, which are on similar 
footing with regard to the string corrections on gravity and electro dynamic sides. We pre- 
sented static topological black hole solutions of third order Lovelock gravity in the presence 
of Born-Infeld gravity, which are asymptotically AdS for negative cosmo logical constant, dS 
for positive A. For the case of solutions with positive curvature horizon (k = 1), one can 
also have asymptotically flat solutions provided A = 0. The topological solutions obtained 
in this paper may be interpreted as black holes with two inner and outer event horizons for 
m > m ext , extreme black holes for m = m cxt or naked singularity otherwise. We found that 
these solutions reduce to the solutions of Einstein-Born- Infeld gravity as the Lovelock coef- 
ficients vanish, and reduce to the solutions of third order Lovelock gravity in the presence 
of Maxwell field as (3 goes to infinity ]6]. We consider thermodynamics of asymptotically 
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flat solutions and found that the first law of thermodynamics is satisfied by the conserved 
and thermodynamic quantities of the black hole. We also consider the rotating solution 
with flat horizon and computed the action and conserved quantities of it through the use 
of counterterm method. We found that the entropy obeys the area law for black branes 
with flat horizon. We obtained a Smarr-type formula for the mass of the black brane as a 
function of the entropy, the charge, and the angular momenta, and found that the conserved 
and thermodynamics quantities satisfy the first law of thermodynamics. We also studied 
the phase behavior of the (n + l)-dimensional rotating black branes in third order Lovelock 
gravity and showed that there is no Hawking-Page phase transition in spite of the angular 
momenta of the branes and the presence of a nonlinear electromagnetic field. Indeed, we 
calculated the heat capacity and the determinant of the Hessian matrix of the mass with 
respect to S, J, and Q of the black branes and found that they are positive for all the phase 
space, which means that the brane is locally stable for all the allowed values of the metric 
parameters. 
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